Multiparticle entanglement under the influence of decoherence 
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We present a method to determine the decay of muhiparticle quantum correlations as quantified by 
the geometric measure of entanglement under the influence of decoherence. With this, we compare 
the robustness of entanglement in GHZ-, cluster-, W- and Dicke states of four qubits and show 
that the Dicke state is most robust. Finally, we determine the geometric measure analytically for 
decaying GHZ and cluster states of an arbitrary number of qubits. 

PACS numbers: 03.67.Mn, 03.65.Yz, 03.67.Pp 



Introduction — The decoherence of quantum states is 
a process in quantum dynamics, which is relevant for the 
discussion of fundamental issues like the transition from 
quantum to classical physics Also from a practical 
point of view decoherence phenomena have to be stud- 
ied, as they occur in experiments involving entanglement 
and their suppression is of vital importance for any im- 
plementation of quantum information processing. 

Due to this importance, the influence of decoherence on 
the entanglement of multiparticle systems has been stud- 
ied from several perspectives 0, Q • These investigations 
concerned either the lifetime of entanglement or the en- 
tanglement properties of the bipartite system which arise 
if the multiparticle system is split into two parts. The 
lifetime of entanglement, however, gives no quantitative 
information about the decay of entanglement More- 
over, as a highly entangled multiparticle state may be 
separable with respect to each bipartition [1], consider- 
ing bipartite aspects only may not lead to a full under- 
standing of the decoherence process. It is therefore highly 
desirable to investigate a full multipartite entanglement 
measure under the influence of decoherence. Unfortu- 
nately, all known entanglement measures for multipar- 
ticle entanglement are defined via complicated optimiza- 
tion procedures , which makes it practically impossible 
to compute them for a given mixed quantum state. 

In this paper we present a method to investigate the 
decay of quantum correlations which can be used to over- 
come these difficulties. We study different four-qubit 
states and use our method to compare their robustness 
against decoherence, using a phenomenological model de- 
scribed below. Our approach allows to compute the en- 
tanglement for GHZ and cluster states of an arbitrary 
number of qubits and thereby to investigate the scaling 
behavior for these states under decoherence. As we will 
further see, our results can be directly tested in nowa- 
days experiments with photons or trapped ions. Finally, 
from the viewpoint of pure quantum information theory, 
our results represent one of the few cases where the com- 
putation of a relevant entanglement measure for mixed 



states can be performed [6|. 

The Situation — We consider the following situation: 
a pure quantum state \^p) is prepared at time t = 0, and 
in the presence of noise evolves to a mixed state g{t). Our 
task is to quantitatively investigate the time evolution of 
the entanglement E{t) = E[g{t)], and its dependence on 
the initial state and the number of qubits. 

As entanglement quantifier we use the geometric mea- 
sure of entanglement This is a popular entangle- 
ment monotone for multiparticle systems, which is re- 
lated to the discrimination of multiparticle states with 



local means 
perspectives 



and has been investigated from several 
MEI. For pure states, it is defined as 



Eg(M) = 1 ^ max 

\'P) = \a)\b)\c). 
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i.e. as one minus the maximal overlap of with fully 
separable states \(f>). It is extended to mixed states by 
the convex roof construction 
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where the minimization is taken over all convex decom- 
positions of g, i.e. over all probabilities pk and states \<pk) 
which fulfill J2kPk\4'k){4>k\ = g- Clearly, the optimization 
in Eq. ^ and especially in Eq. ^ is difficult to perform. 

Our method can be summarized as follows: since any 
set of probabilities p/c and states in Eq. ^ results in 
a valid upper bound, we obtain a good upper bound by 
choosing them appropriately. Then we use the results of 
Refs. [12,ll3l to obtain a lower bound on Ea(t). There it 



has been shown how the geometric measure can be esti- 
mated if the mean value of a single or a few observables is 
given 1J| . We show that the lower and the upper bound 



coincide for the multiparticle states we investigate below, 
allowing for a precise determination of Ecit). 

The noise we consider is described by a master equa- 
tion for the matrix elements gki as they are used in 
phenomenological models of decoherence for e.g. electron 
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spin qubits 



dtQki = 
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WkiQkk) 



for k ~ I, 
for k^l. 



(3) 



We consider a global dephasing process, where the relax- 
ation of the diagonal elements plays no role {Wij = 
Vi, j) and the off-diagonal terms are affected by a global 
dephasing rate Vki = 7. This process leads to exponen- 
tially decaying off-diagonal components gki{t) = xgki{0), 
where here and in the following x = e"'''*. At the end 
of the paper we will discuss extensions of our method to 
other (e.g. local) decoherence models. 

Four-qubit states — Before presenting our results for 
different four-qubit states in detail, note that by con- 
struction the geometric measure is a convex quantity, 
i.e. EaiXgi + (f - X)g2] < XEaigi) + (1 - \)EG{g2). 
From this, it follows directly that Ecit) is monotonically 
decreasing. Moreover Eq{x) is convex in the parameter 
X [since g{x) depends linearly on a;], and consequently 
EQ{t) is convex in the time t. 

Let us start our discussion with the four-qubit 
Greenberger-Horne-Zeilinger (GHZ) state [l^, given by 



\GHZi) = (10000) + |lflf 



(4) 



The geometric measure of the GHZ states equals 1/2 as 
the maximal overlap with product states is 1/2 

In order to estimate the entanglement from below, 
note that the time dependent fidelity is given by Fit) = 
Tr[g{t)\GHZ4){GHZ4\] = (l-he-''*)/2. Using the results 
of Ref. we can estimate the geometric measure from 
F{t). Explicitly, it has been shown that if for an arbi- 
trary g the fidelity F of a state 1-00) is given, then Eq is 
bounded by 



Ecig) > sup{l+rF(t)-i[l+r+v/(l + r)2 - irEai^o)]}- 

r>0 

(5) 

Applying this to g{t) and the fidelity of the GHZ state 
leads to (using x — e^^*) 



Ecig) > ^i-Vi-x^)- 



(6) 



In order to obtain an upper bound, we consider the 
two states |0i) = c|0000) + s|llll) and |</)2) = s|0000) -I- 
c|llll) with c = cos(a),s = sin(a), and write g{t) = 
(l/2)Efc=i,2 \<l^k){4>k\- Then, x/2 = (2cs)/2 has to hold 
and using the fact that the geometric measure for the 
states \(j}k) is given by Eq = minjs^, c^} [13] one arrives 
at an upper bound for Ecit) which is given by the right 
hand side of Eq. ^. Hence, Eaig) = (1 - \/l~x^)/2 
for the four-qubit GHZ state ^ under the influence of 
noise. This function is shown in Fig. [TJ 

Second, let us discuss the four-qubit cluster state [l^, 

\CLi) = (10000) + 10011) + |1100) - |llll))/2, (7) 
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FIG. 1: The geometric measure Ecit) for four-qubit GHZ-, 
cluster-, W- and Dicke states for the case 7 = 1. For the W- 
and Dicke state, the upper bounds are shown, as the curves 
of the lower bounds coincide with that. For t > 2 the values 
for the W- and cluster state coincide [see Eqs. (f8l llOp]. 

which has a geometric measure of i?G(|C-^4)) = 3/4 [lo| . 
For this state, the decay of the fidelity is given by F{t) = 
(1 + 3a;)/4, which [combined with Eq. ^] leads to the 
lower bound 



Eig) > I [1 + X - y/l + {2- 3x)x] . 



(8) 



For the upper bound, we consider four trial vectors 
for the decomposition. The first is given by = 
c|0000) -I- s(|1100) 4- [0011) - |llll))/\/3 and the other 
three are obtained from this by permuting the four terms. 
We choose c > s, then any of the four states has a 
geometric measure of EG{\4>i)) = s^. With the ansatz 
g{t) — (l/4)X]fc=i \4'k){4>k\ we obtain as condition on s 
and c that x/A = [{2cs/V3) + (2s2/3)]/4. From this, c 
can be determined. This leads after a short calculation 
to the insight that the right hand side of Eq. ([8]) also con- 
stitutes an upper bound on the entanglement, and thus 
describes exactly the time evolution of the entanglement. 
Third, a four-qubit W state is given by [T^ 

IW4) = (10001) + [0010) + [0100) + |1000))/2. (9) 

for which the geometric measure equals 37/64 Let us 
first derive the upper bound. We take as test states the 
state |0i) = c|1000) + s(|0100) + [0010) + |0001))/V3 
and permutations thereof. Using a symmetry argu- 
ment 0, their geometric measure is determined to be 
Ecim) = [5 + 3c2(c4 + c2-3)]/[(3-4c2)2l for c < 
I/V2 and EG{\(t)^)) = 1 - c2 for c > 1/^2 [l^- Then, 
one can derive an upper bound as for the cluster state. 

It turns out, however, that this upper bound is not 
convex in x. Since we know that Eg has to be convex, 
we can take the convex hull (in x) of this upper bound, 

37(81a;-37) , 
Eg < ^ ^g^g ' for X > xo, (10) 

Eg < I [1 + 2^ - Vl + (2 - 3a;)a;] for x < xq, 



3 



with Xq = 2183/2667. Physically, taking the convex hull 
just means that for short times (when x > xq) the opti- 
mum in the convex roof in Eq. ^ is of the form g{x) = 
p|W^4)(VF4| + (l-p)^?(xo),with gixo) = {l/4)Ek\^k){^k\- 
Note that for longer times {x < xq) the upper bound HU]) 
is the same as for the cluster state. 

In order to see that this upper bound is optimal, let 
us derive a lower bound. Here, we not only take the fi- 
delity of the W state into account, but we use as a second 
observable the projector onto the space with one excita- 
tion, Vi = 10001) (0001| -I- 10010) (00101 + 10100) (01001 -I- 
|1000)(1000|. Using the fact that the fidelity of IW4) is 
given by F{t) = (1 -I- 3a;)/4 and that it is always in the 
space spanned by (i.e., Tr[g{t)Vi] — 1), we can use 
the methods of Ref. [ij] to obtain a lower bound from 
these two expectation values [2^. It turns out that, 
within numerical accuracy, the lower bound coincides 
with the upper bound, giving strong numerical evidence 
that Eq. (fTO|) is the exact expression for the decay of 
quantum correlations in the W state ([9]). 

As a last example of a four-qubit state let us discuss 
the symmetric Dicke state 2l| given by 

\D4) = (10011) + 10101) + |1001) + 

+ 11100) + 10110) -I- 11010))/\/6, (11) 

which has a geometric measure Eq = 5/8 Q- To obtain 
the upper bound, we proceed similarly as for the W state: 
We take six states, the first one being = c|0011) -I- 
s(|0101) |1001) + 11100) + 10110) + \W10))/V5 and other 
ones obtained by permuting the terms. Then we make 
the ansatz g{t) — {^/6)J2k=i \4'k){(f'k\- This leads to an 
upper bound which is not convex in x, and subsequently 
taking the convex hull leads to 



5(3x-l) ^ 5 
Eg < — for a; > -, 



(12) 



Eg < ^ [1 -t- 22; - Vl + (4-5x)x] for a; < f . 

The lower bound is found analogously to the W state 
as well, with the projector V2 onto the space with two 
excitations as second observable. The resulting bound 
coincides again with the upper bound, giving strong ev- 
idence that Eq. ([T^ describes the time evolution of the 
entanglement. 

Comparison of the four-qubit states — For this com- 
parison we consider the logarithmic derivative ry = 
i9t(ln[i?G(i)]) = dtEG(t)/EG(t) which describes the rela- 
tive decay of entanglement [2^ . The values of this quan- 
tity are plotted in Fig. [21 One can clearly see, that the 
Dicke state is the most robust state, while the GHZ state 
is the most fragile state [2^. It is an interesting open 
question which properties of the Dicke state are respon- 
sible for the high robustness. 

Many qubits — We restrict our attention to GHZ and 
linear cluster states, as they are highly relevant for ap- 
plications like quantum metrology or measurement-based 
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FIG. 2: Logarithmic derivative of Eait) (for 7 = 1) for dif- 
ferent four-qubit states. The non-analytic behavior for the 
Dicke- and W state stems from the convex hull in Eqs. (|10I12|I . 



quantum computation 2J|. In the decoherence model, 
one might keep the dephasing rate 7 = 70 constant for 
any number of qubits (where 70 is the single qubit de- 
phasing rate), or scale it as 7 = iV7o (as it would oc- 
cur for the GHZ state in a local dephasing model). For 
the GHZ state, nothing changes and all the formulae ob- 
tained in the above for the four-qubit case apply. Con- 
cerning the cluster state, we consider linear cluster states 
with N = 2n qubits. The linear cluster state is given by 



\CL 



N 



+ \n){a,(g>l)]/V2, (13) 



fc=i 



where this formula should be understood as an itera- 
tion, with the operator {a^ <8) 1) acting on the Bell state 
of the next two qubits. Explicitly, we have ICL2) = 
(|00)-f |ll))/\/2and \CLi) - |00)(|00)-f |ll))-H|ll)[a,0 
1(|00) + 111))] = 10000) + 10011) + |1110) + |1101) 0- 
Note that the maximal overlap of the cluster state with 
fully separable states is 1/2" [10]. 

We calculate the geometric measure for this state un- 
der the effect of decoherence in the same way as for the 
four-qubit case. The lower bound is obtained from the fi- 



delity F(i) = [l + (2'^ 
3 • 2" -I- 2^)x + 2(2" ■ 



~l)x]/2" and yields £^0 > 2^((2- 
l){l - y(r^^)[rTT2^^^^1>I}). 
For the upper bound, we consider 2" test states similar 
to the ones before, and arrive at an upper bound which 
coincides again with the lower bound. 

To investigate the scaling behavior, we consider the 
time ti/2, when the entanglement has decreased to half 
of the initial value. These times can be directly computed 
for the A^-particle GHZ and cluster states. Fig. [3] shows 
ti/2 as a function of the number of qubits N, for the con- 
stant and linear models of 7. In both cases the time ti/2 
of the cluster state is, in the limit N 00, larger than 
that of the GHZ state by a factor of ln(4)/ ln(4/3) « 4.82, 
giving quantitative evidence for the higher robustness 
against dephasing of the cluster state. 
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FIG. 3: Comparison between the GHZ state and the cluster 
state as a function of the number of qubits A*'. The half life 
time ti/2 is shown for the case that 7 = 4 and for the case 
that 'y = N increases linearly with the number of qubits. 



Discussion — In the calculations presented in this 
paper we concentrated on a global decoherence model. 
However, our results can also be applied to other models. 
First, for the W- and GHZ state, our model is equivalent 
to a local dephasing noise, as it occurs in multi-photon 
experiments. Given the experimental availability of W- 
and GHZ states with photons j20| , our results can be di- 
rectly tested with present-day technology. Second, for lo- 
cal decoherence models where relaxation is the dominant 
process, a small modification of our scheme allows the 
calculation of the geometric measure for certain states, 
such as the W state 27| . The fact that this type of deco- 
herence is dominant in ion traps 12 8| combined with the 
possibility to generate such states [29| opens another way 



for an experimental test. 

To conclude, our results provide a novel and versatile 
method to determine the decay of multiparticle entan- 
glement for quantum states under the influence of deco- 
herence. Our results can be tested in multi-qubit experi- 
ments and may therefore lead to a better understanding 
of decoherence as a fundamental obstacle in quantum in- 
formation processing. 
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